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We report thermopower measurements of interacting GaAs bilayer hole systems. When the carrier 
densities in the two layers are equal, these systems exhibit a reentrant insulating phase near the 
quantum Hall state at total filling factor v = 1. Our data show that as the temperature is decreased, 
the thermopower diverges in the insulating phase. This behavior indicates the opening of an energy 
gap at low temperature, consistent with the formation of a pinned Wigner solid. We extract an 
energy gap and a Wigner solid melting phase diagram. 
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The ground state of two-dimensional (2D) carrier sys- 
tems has been subject to intensive work over the past 
decades. The interplay of Coulomb interaction, disor- 
der, and magnetic field B leads to many different ground 
states 0. In low-disorder 2D systems, a sufficiently large 
B is expected to quench the kinetic energy of the carri- 
ers, leading to the formation of a Wigner solid (WS) 0. 
Experimental evidence has been accumulated for the 
existence of a WS in both GaAs 2D electron single- 
layers 0, 0, around the fractional quantum Hall 
state at Landau level filling factor v = ^ and electron 
bilayers in the vicinity of total filling factor v — i , as 
well as in dilute, single-layer, GaAs 2D hole systems |fj 
near v = i. 

Recently, a reentrant insulating phase (IP) was ob- 
served in interacting GaAs bilayer hole systems around 
the many-body quantum Hall state (QHS) at total fill- 
ing factor v = 1 0. It was suggested in Ref. that 
this IP signals a bilayer WS, stabilized at such high v by 
strong Landau level mixing due to higher effective mass 
of holes, and by the interlayer Coulomb interaction in 
the bilayer system. While magnetotransport studies in 
bilayer systems have intensified in recent years 0, lit- 
tle attention was paid to their thermoelectric properties. 
Thermopower coefficients are sensitive tools to study the 
electronic properties of 2D systems, including Si inversion 
layers 0, SiGe quan tum wells 10], and GaAs hetero- 
junctions [111 Il2t In particular, as the temperature 
T — > 0, the diffusion thermopower is expected to van- 
ish as T 1 / 3 for a mobility gap or diverge as T _1 in the 
presence of an energy gap |l4| . 

We study here the thermopower of interacting GaAs 
bilayer hole systems around v = 1. When the hole 
densities in the two layers are equal, at filling factors 
near v = 1, our measurements show that the diffusion 
thermopower diverges as T™ 1 when the temperature is 
lowered below ~ 100 mK. This behavior corroborates 
the magnetotransport and the non-linear / — V re- 
sults. The T dependence of the thermopower indicates 
the opening of an energy gap, consistent with the forma- 



tion of a pinned, bilayer WS near the QHS at v = 1, and 
allows us to extract a phase diagram for the WS melting. 

The samples we have studied are Si modulation-doped 
GaAs bilayer hole systems grown by molecular beam epi- 
taxy on undoped GaAs (311)A wafers. The two samples 
consist of two 15 nm wide GaAs quantum wells separated 
by an AlAs barrier with a thickness of 11 nm (sample A) 
and 7.5 nm (sample B). The width of the AlAs barrier 
and the large effective mass of the holes assure a weak in- 
terlayer tunneling [To] ] . Two 8x2 mm 2 rectangular pieces 
of wafer were cut with their length along the [011] (sam- 
ple A) and [233] (sample B) crystallographic directions. 
On both specimens, we patterned 400 /jm wide Hall bars 
in order to reduce density inhomogeneities; the distance 
between the voltage probes along the Hall bar is 1.2 mm. 
Ohmic contacts were made to both layers. A back-gate 
was used to tune the charge distribution of the two layers. 
The charge imbalance is denoted by Ap = (px — Pb)/2, 
where ps and pr are the hole densities of the bottom and 
the top layer, respectively. For balanced charge distribu- 
tions (Ap = 0), the total hole density for both samples is 
close to p — 5 x 10 11 cm -2 [J^. Our experimental s etup 
for thermopower measurements is described in Ref. |12j . 
One end of each sample was mounted onto the cold fin- 
ger of a dilution refrigerator, while a strain gauge, used 
as heater, was glued on the other end using GE varnish. 
To measure the longitudinal thermopower S xx , a ther- 
mal gradient was created along the sample by applying a 
sine wave current at a frequency / = 1.3 Hz through the 
heater and measuring, via a lock-in amplifier, the voltage 
signal at 2/ between two contacts on the same edge of 
the Hall bar. The temperature gradient was determined 
using two carbon paint thermometers deposited on the 
sample's surface, and was kept below 10% of the mean 
T. 

We first summarize the magnetotransport results for 
our samples. In Fig. ^ we illustrate for sample A the 
longitudinal resistivity p xx as a function of B at low T. 
Consistent with the sample's extremely small interlayer 
tunneling 
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at low B the p xx trace exhibits minima 
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FIG. 1: Upper panel: longitudinal resistivity p xx as a function 
of B at the indicated temperatures for sample A at a total 
density p = 5.1 x 10 11 cm -2 . For the data in the figure, the 
densities of the individual layers are equal (Ap = 0). Top 
inset: The T dependence of p xx at v = 0.92, 0.95, and 1.1. 
Bottom inset: A magnified view of the low-B p xx trace at 
30 mK. Lower panels: differential resistance in the insulating 
phase of sample A as a function of the dc bias for v = 1.1 
(left panel) and 0.95 (right panel) at 30 mK. 



only at even integer fillings v = 2, 4, 6, etc. The strong 
interlayer interaction, however, leads to a many-body 
QHS at v = 1 ||, flanked by reentrant IP's 7]. The 
measured QHS excitation gaps : A, extracted from the 
T dependence of p xx at v — 1, are 0.13 and 0.67 K, re- 
spectively; the larger 1 A for sample B is consistent with 
its smaller d/ls ratio 0|. The T dependence of p xx at 
v = 0.92, 0.95, and 1.1 is displayed in the top inset to 
Fig. n revealing the insulating behavior at these fillings. 
Finally, in the bottom panels of Fig. ^ we show I — V 
characteristics at filling factors near v = 1 for sample 
A, exhibiting a strongly non-linear behavior. We have 
observed qualitatively similar I — V curves for sample B 
also. These non-linear data are similar to previous mea- 
surements in the IP near v — | in single-layer GaAs 2D 
hole systems 

Before discussing our S xx data, we briefly review the 
essential concepts for understanding thermopower data 



in 2D carrier systems . In general S xx consists of two 
additive contributions: the phonon-drag thermopower 
S!j. x and the diffusion thermopower S xx . In GaAs-based 
hetcrostructures, the former is known to be dominant for 
T > 200 mK, while the latter provides the main contribu- 
tion at lower temperatures (T < 100 mK). The phonon- 
drag thermopower is given by S% x — ^: v ^. h , where C 
is the specific heat, A is the phonon mean free path, e 
is the electron charge, v s is the sound velocity, and T p h 
is the phonon momentum relaxation time due to scat- 
tering with holes only [l7j] . The diffusion thermopower, 
on the other hand, measures the heat or the entropy per 
carrier and is sensitive to the density of states at the 
Fermi level. In particular, for a 2D insulating system, 
S xx can distinguish between a mobility gap and an en- 
ergy gap. In the first case, S xx should vanish at low T as 
S xx oc T 1 / 3 , while for the second case - when an energy 
gap Eg separates the ground state from its excitations - 
S xx is expected to diverge at low T as 
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Hence, S xx is a unique probe of the ground state of 2D 
systems as it can critically test if the IP results from the 
strong localization of the carriers (mobility gap) or from 
the formation of a pinned WS (energy gap) . 
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FIG. 2: S mx as a function of B at four temperatures below 
120 mK for sample A. Inset: S X x vs B is displayed at T = 285, 
310, and 365 mK in the phonon-drag regime. 

We now present the result of our S xx measurements, 
summarized in Figs. |3 and and consider first the 
phonon-drag regime. Figure [5] shows S xx vs B traces 
for sample A at various temperatures. The curves at 
T = 285, 310, and 365 mK, where the phonon-drag is 
dominant, are displayed in the Fig.|21inset. Similar to p xx 
data, we observe minima at even filling factors (v = 2, 
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4, and 6) and at v = 1 consistent with a vanishing den- 
sity of states at the Fermi level ^3 ■ -"- n the phonon-drag 
regime, a dramatic decrease of S xx is found in the en- 
tire investigated B range as T diminishes. This is clearly 
indicated by the T evolution of S xx at various fillings 
shown in Fig. [3] for both samples. In the top panel of 
Fig. S xx vs T is displayed for sample A at v = 1.1, 
1, 0.95, 0.92 and for B = 0.1 T when the bilayer system 
is balanced. In the bottom panel of the figure, we show 
S xx vs T for sample B. Data points are presented at v — 
1.14 for a balanced and a slightly imbalanced charge dis- 
tribution. At high T, S xx = S!^ x follows the well known 
T 4 law 201 for all the experimental configurations shown 
in Fig. 111. 
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FIG. 3: Top panel (sample A): T dependence of S xx at differ- 
ent filling factors when the bilayer is balanced. Bottom panel 
(sample B): T dependence of S xx at u = 1.14 for balanced (•) 
and imbalanced (o) charge distributions. The dotted lines in- 
dicate various T dependencies. The error bars are given by 
the noise of the lock-in amplifier, the magnitude of the applied 
thermal gradient, and the accuracy of the thermometers' cal- 
ibration. 



in Fig. Focusing on the data shown in the top panel of 
Fig. we observe that at v = 1, S xx rapidly ap pro aches 
zero as expected for a QHS with no entropy [lj] and 
S xx at B = 0.1 T is proportional to T below 200 mK 
as expected for a diffusion-dominated thermopower in a 
metallic system. However at v = 1.1, 0.95, and 0.92, 
S xx oc T _1 as T — > 0. This divergence of S xx in the 
reentrant IP's at very low T indicates the opening of an 
energy gap and suggests the formation of a pinned bilayer 
WS near v = 1. 

Data for sample B exhibit qualitatively the same be- 
havior. In the lower panel of Fig. [31 we show S xx vs T 
data for sample B at v = 1.14 where this sample's mag- 
netotransport data reveal an IP. The S xx data show a 
divergence at low T when the bilayer system is balanced 
(closed circles, Ap = 0). In the same panel we also show 
data (open circles, Ap = 2.6 x 10 9 cm -2 ) for sample B 
when charge distribution is imbalanced, i.e., when charge 
is removed from the the bottom layer via the application 
of positive back-gate bias. For the imbalanced state, S xx 
at v = 1.14 turns from a T~ x to an approximatively 
linear temperature dependence. This observation is in 
agreement with magnetotransport results showing that 
the reentrant IP disappears for small charge imbalance 
between the two layers Q. We note that the strength- 
ening of the v — 1 QHS when Ap increases, seen in the 
magnetotransport is also observed in thermopower 
measurements. In other words, the filling factor range 
around v = 1 where S xx — > as T — > extends when 
charge is transfered from one layer to the other (data not 
shown here). 
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At low temperatures, where the thermopower is dom- 
inated by the diffusion contribution, S xx increases with 
decreasing T in the IP's reentrant around v = 1. This 
main result of our work is illustrated in Fig. El by the S xx 
vs B traces at T = 113, 71, 58, and 51 mK 22]. The 
divergence of S xx at low T in the IP is further illustrated 



FIG. 4: Evolution of the melting temperature T m as a func- 
tion of v of the pinned, bilayer WS in sample A. The top 
axis shows the corresponding magnetic fields. Inset: T m is 
extracted from the T dependence of S xx by assuming that 
the rise in S xx signals the melting of the WS. 
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We now turn to a more quantitative discussion of our 
thermopower data. From the T dependence of S xx , 
shown in Fig. [3] and the inset to Fig. 21 we can extract a 
phase diagram for the WS melting. To estimate the melt- 
ing temperature T m , we assume that the solid-to-liquid 
transition occurs when S xx starts to rise as T is low- 
ered. As illustrated in the inset to Fig. 0] for v = 0.93, 
we have fitted at fixed filling factors the phonon-drag 
thermopower with a T 4 law for T > 200 mK and the 
diffusion thermopower with a T" 1 law for T < 100 mK. 
The intersection of these two fits determines T m (y). In 
Fig. 01 the T m vs v phase diagram is shown for sample 
A. The reentrant behavior of the IP is illustrated by the 
evolution of T m ; T m is finite in a small filling range cen- 
tered near v = 1.1 and at v > 0.95, but is vanishingly 
small near v = 1 (B = 2.1 T) where the quantum Hall 
ground state is present. The magnitude of the estimated 
T m in our sample is similar to the reported values of 
the WS melting temperature in single-layer 2D hole 0] 
and electron pi systems around v — \ and |, re- 
spectively. Using Eq. and the S xx = aT^ 1 fit to the 
diffusion regime, an estimate of the energy gap associated 
with the pinned WS can be made when E g 3> k B T [l2| . 
This energy gap likely corresponds to the creation en- 
ergy of thermally activated defects (such as bound pairs 
of dislocations), which may be responsible for electri- 
cal conduction in the pinned WS |2J|. For example, at 
v = 0.92, E g /k B = 0.34 K in sample A. This value is 
of the same order of magnitude as the activation energy 
^a/^b = 0.15 K determined from fitting p xx in the low- 
T range to ln(p xx ) oc E^/2k B T, and compares well with 
the WS melting temperature T m = 0.18 K at this filling. 

To conclude, we have measured the thermopower of 
interacting GaAs bilayer hole systems around v = 1. At 
very low T, when the diffusion thermopower dominates, 
S xx diverges in the IP's reentrant around v = 1. This 
behavior indicates the opening of an energy gap and the 
presence of a pinned bilayer WS. We have interpreted 
our thermopower data in terms of a WS melting and 
extracted a T m vs v phase diagram. 
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